System and method for retrieving an image containing image statistical data

ABSTRACT

A system and method of retrieving an area of interest of an image stored in a storage device. The method includes selecting the area of interest having user desired features, the stored image being divided into a plurality of polygonal image sections, each polygonal image section having substantially homogenous features, the image being further divided into a plurality of quads such that each quad has boundaries that contains at least a portion of a first polygonal image section and at least a portion of a second polygonal image section. The method further includes determining in which quad in the plurality of quads each pixel of the area of interest is located and determining for each quad, in which portion of the first polygonal image section or the second polygonal image section each pixel is located.

CROSS REFERENCE TO RELATED APPLICATIONS

The present patent application is a continuation application of U.S. patent application Ser. No. 13/779,251 files on Feb. 27, 2013 which is a continuation of U.S. Pat. No. 8,411,970 issued on Apr. 2, 2013 which claims priority from provisional patent application No. 61/314,257 filed on Mar. 16, 2010, the entire content of each are incorporated herein by reference.

BACKGROUND OF THE INVENTION

Field of the Invention

The present invention pertains to data management and in particular to a method and system for determining statistical data for pixels of an image and a method and system of converting an image for viewing on a display device.

Discussion of Related Art

An image generally contains a plurality of pixels and each pixel has one, two or more bands. Each band has a certain color depth or bit depth. For example, an RGB color-based image has 3 bands, the red band (R), the green band (G) and the blue band (B). Each of the R, G and B bands can have a depth of 8 bits or more. In order to be able to visualize or view images having a bit-depth greater than 8 bits per band (i.e., images having a depth of N-bits per band, with N greater than 8) on a display device that is configured to display an image having 3 bands (each band having an 8-bit depth), a bit depth down-sampling operation is performed. For example, a bit-depth down-sampling operation is performed on the image having N-bits per band to transform or convert the image into another image in the RGB system where each band (R band, G band and B band) has 8 bits of depth (i.e., a 24-bit RGB image).

Some conventional methods use a pre-computed N-bit to 8-bit look up table (LUT) per band to down-convert the image with higher bit depth per band into a lower bit depth (e.g., 8-bit depth) per band. When referring to “higher” bit depth per band herein, this is intended to refer to a bit depth of greater than 8 bits per band. Other methods such as tone mapping method involve executing algorithms to adaptively compute the N-bit (where N is greater than 8) to 8-bit conversion, pixel by pixel. The latter method of performing the down-conversion or down-sampling can be computer and time intensive based on the quality of the implementation of the algorithm. An increasing number of computing environments have implemented these two types of down-conversion methods into hardware or into software. Most of these methods employ the generation and application of statistical information gathered by processing pixel data in the image.

When pixels containing bands with bit depth greater than 8-bits per band are stored into disk files, they are stored in 16-bit, 32-bit, or higher bit depths per band. For example, an original image sensor may have captured data at 10, 11, 12, or 15 bits per band. However, when an image file containing these pixel values is stored on disk, it is generally stored at 16-bits per band. While this may be considered as waste of space, it is also deemed highly efficient for computing platforms of today.

When a vendor distributes to a customer imagery that was originally captured at a value greater than 8 bits per band (e.g., 16-bits per band), very often the vendor is forced to distribute the original image data as well as additional image data that was pre-converted at significant computational cost to the 8-bit depth per band, specifically for display. Under these circumstances, the 8-bit data would occupy an additional 50% more storage alongside the original 16-bit data. Compression is further used to save space on disk. If the compression is lossy, the compression degrades the image. If the compression is not lossy, the compression is not as efficient to read and decode the compressed data for visualization. A customer sometimes may choose to discard the original 16-bit data and just use the 8-bit data in critical decision making tasks because of a lack of sufficient amount of storage or the lack of computing power to handle on-the-fly conversions. If the customer chooses to only use 16-bit data in critical decision making tasks, the overhead of on-the-fly conversion to 8-bit depth for visualization bears an additional computational cost or requires additional computing power which further increases costs. However, the flexibility that on-the-fly 16-bit to 8-bit depth conversion offers is that the customer can select one of several methods to convert the image to 8-bit depth prior to display.

If a computationally inexpensive method was available for better on-the-fly visualization of data that was collected and stored at bit depths greater than 8-bits per band, it would be ideal to supply the original data (greater than 8-bits per band) to the customer and let the customer choose amongst the various methods available for displaying this data. Hence, distributing pre-converted 8-bit depth per band data would not be necessary, thus saving additional disk space requirement as well as saving the vendor the added cost of converting the data to 8-bit depth.

The method or methods described in the following paragraphs provide a mechanism for pre-computing and storing data within the image, and utilizing this data to make subsequent on-the-fly conversion of the image data having a higher bit depth per band to a lower bit depth per band (e.g., 8 bit depth per band) for the display of the image computationally inexpensively.

BRIEF SUMMARY OF THE INVENTION

An aspect of the present invention is to provide a method for determining statistical data for pixels of an image, the pixels having a higher bit depth per band. The higher bit depth per band includes greater than 8 bits per band. The method includes dividing the image having pixels with the higher bit depth per band into a plurality of polygonal sections, each polygonal section including substantially homogenous features; computing statistical data for pixels within each polygonal section; and storing the statistical data with the image.

Another aspect of the present invention is to provide a system for determining statistical data for pixels of an image, the pixels having a higher bit depth per band. The higher bit depth per band includes greater than 8 bits per band. The system includes a processor configured to: (a) divide the image having pixels with the higher bit depth per band into a plurality of polygonal sections, each polygonal section including substantially homogenous features; (b) compute statistical data for pixels within each polygonal section; and (c) store the statistical data with the image.

Another aspect of the present invention is to provide a method of converting a first image within an image having pixels with a relatively higher bit depth per band into a second image having pixels with a relatively lower bit depth per band. The method includes selecting a first image, the first image corresponding to a viewport from an image having pixels with the higher bit depth per band, the image having pixels with the higher bit depth per band being divided into a plurality polygonal sections, each polygonal section including substantially homogenous features; and determining, for each pixel with the higher bit depth per band within the first image, to which polygonal section in the plurality of polygonal sections each pixel belongs to. The method further includes converting each pixel with the higher bit depth per band within the first image into a corresponding pixel of a lower bit depth per band using statistical data of a polygonal section to which each pixel belongs to; and forming the second image having pixels of the lower bit depth per band.

Yet another aspect of the present invention is to provide a system of converting a first image having a relatively higher bit depth per band into a second image having a relatively lower bit depth per band. The system includes a processor configured to: (a) select a first image, the first image corresponding to a viewport from an image having pixels with the higher bit depth per band, the image having pixels with the higher bit depth per band being divided into a plurality polygonal sections, each polygonal section including substantially homogenous features; (b) determine, for each pixel with the higher bit depth per band within the first image, to which polygonal section in the plurality of polygonal sections each pixel belongs to; (c) convert each pixel with the higher bit depth per band within the first image into a corresponding pixel of a lower bit depth per band using statistical data of a polygonal section to which each pixel belongs to; and (d) form the second image having pixels of the lower bit depth per band.

Although the various steps of the method are described in the above paragraphs as occurring in a certain order, the present application is not bound by the order in which the various steps occur. In fact, in alternative embodiments, the various steps can be executed in an order different from the order described above or otherwise herein.

These and other objects, features, and characteristics of the present invention, as well as the methods of operation and functions of the related elements of structure and the combination of parts and economies of manufacture, will become more apparent upon consideration of the following description and the appended claims with reference to the accompanying drawings, all of which form a part of this specification, wherein like reference numerals designate corresponding parts in the various figures. In one embodiment of the invention, the structural components illustrated herein are drawn to scale. It is to be expressly understood, however, that the drawings are for the purpose of illustration and description only and are not intended as a definition of the limits of the invention. As used in the specification and in the claims, the singular form of “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise.

BRIEF DESCRIPTION OF THE DRAWINGS

In the accompanying drawings:

FIG. 1 depicts the geographical map of the continental United States of America (CONUS) that is divided into visually homogenous regions or sections, according to an embodiment of the present invention;

FIG. 2A depicts a converted image from an unsigned 16 bits per band image using a conventional standard-deviation technique applied to the histogram of the whole image file;

FIG. 2B depicts the same image converted from unsigned 16 bits per band to an unsigned 8 bits per band image according to an embodiment of the present invention;

FIG. 2C depicts a polygonal layout (e.g., quadrilateral layout) used to sub-divide the image in FIG. 2B into a plurality of regions or sections; according to an embodiment of the present invention;

FIG. 3A depicts an image portion obtained when zooming into the top-left corner region of the image shown in FIG. 2A;

FIG. 3B depicts an image portion obtained when zooming into the top-left corner region of the image shown in FIG. 2B;

FIG. 4 is a representation of R-sets in relation to each other and in relation to the original image, according to an embodiment of the present invention;

FIG. 5 represents a hierarchical tree or relationship structure between various quads within an image, according to an embodiment of the present invention;

FIG. 6 depicts a structure for storing statistical data for the original image and each R-set, each divided into several regions, according to an embodiment of the present invention;

FIG. 7 is a flow chart depicting a method of pre-computing the statistical data used for converting a first image of higher bit depth to a second image of lower bit depth, according to an embodiment of the present invention; and

FIG. 8 is a flow chart depicting the method of converting a first image to a second image that can be displayed for a given area of interest (AOI) in the NUI file, using feature-based statistics and associated data generated earlier, according to an embodiment of the present invention.

DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION

Some image files such as image files in the new universal image (NUI) file format implemented by PIXIA Corporation, can contain a relatively large number of pixels that can scale to over four billion (4×10⁹) by four billion (4×10⁹) pixel raster. NUI image files having a size of several million pixels wide by several million pixels tall can be found in many areas including geography, aerial and satellite images, medicine, etc. For example, images that cover many geographical areas such as cities, provinces, countries and continents can be as large as several million pixels by several million pixels. The NUI file format has the capacity of storing typical band layouts such as luminance, red-green-blue (RGB), near infrared (NIR), multi-spectral, hyper-spectral and ultra-spectral bands. The NUI file format can also be configured to store various bit-depths such as, but not limited to, unsigned 8 bits (u8), signed 16 bits (s16), unsigned 16 bits (u16), signed 32 bits (s32), unsigned 32 bits (u32), 32-bit floating point (f32), 64-bit double precision floating point (f64) and complex. In addition, within NUI files containing bands with signed 16-bit depth and unsigned 16-bit depth, the NUI file format can also carry between 9 and 15-bits per band image data. The NUI file format can store the image data uncompressed. Alternatively, the NUI file format can also store the image data as compressed. In one embodiment, the compression can be performed using a conventional JPEG2000 compression method.

The NUI file, as part of its format, stores statistical data about the image. This statistical data includes, but is not limited to, a histogram of each band of each pixel in the image, minimum, maximum value of each band of each pixel in the image, standard deviation, mean, median, mode, skew, and kurtosis of each band of each pixel in the image. In one embodiment, a display system or device can use the statistical data to generate a custom look-up table (LUT). Multiple LUTs may be generated using multiple techniques based on user's requirements. The LUT may be further used to convert pixels in the image containing bands of higher bit depth (e.g., bit depth greater than 8-bits per band) to an image containing bands of lower bit-depth (e.g., 8-bits per band) that can be displayed on, for example, a 24-bit RGB display device (i.e., a display device limited to an 8-bit depth per band).

The NUI file, as part of its format, stores in addition to the original image, a set of user-defined reduced resolution datasets (RRDs or R-sets). Each RRD is a smaller, zoomed-out version of a previous RRD and a first RRD is a smaller, zoomed-out version of the original image. The terms RRD and R-set are used herein interchangeably. FIG. 4 is a representation of R-sets in relation to each other and in relation to the original image, according to an embodiment of the present invention. As shown in FIG. 4, the first R-set is a zoomed-out version of the original image, the second R-set is a zoomed-out version of the first R-set, and the third R-set is a zoomed-out version of the second R-set, etc.

However, the use of the histogram for the entire NUI file, in particular when the NUI file is relatively large (e.g., a file containing imagery captured over a vast geographical region such as satellite or aerial imagery), may have a drawback. Indeed, as the number of populated pixels in the NUI file increases, the dynamic range of pixel band color values stored within the image data can increase considerably. For example, a large geographical area occupied by a country such as the United States of America, China, Russian Federation or India may cover various geographical features such as deserts, tropical forests, farmlands, snow capped or rocky mountain ranges, and densely populated cities. Therefore, a histogram obtained from such a diverse set and vast region of pixels would have a very broad distribution of color (e.g., golden, brownish, green, blue, etc.).

If such a diverse dynamic range of pixels is converted from higher bit-depths, for example from unsigned 16 bits (u16 bits) per band to unsigned 8 bits per band within a 3 band color scheme, a single u16 to u8 look-up table (LUT) may not have enough slots to represent all the visual features and the range of colors within these features.

The resulting converted image (converted from higher bit depths per band to lower bit depths per band) may look “washed out” or “faded” with every region losing detail so as to do distribute the limited color space to all other regions. This is generally rectified by applying further image processing methods such as a contrast curve or dithering to the converted image. Since the detail is already lost, the best this process can do is add additional artifacts that further remove detail but give the perception of greater image fidelity. For example, subtle gradations in clouds and snow covered regions of very rocky mountains may not show up and would look to the eye all white. This may not be desirable when implementing high dynamic range rendering (HDRR) and tone mapping to retain a high dynamic range image (HDRI). Hence, in order to remedy this deficiency, in one embodiment, a plurality of histograms can be generated and stored for geographical regions that contain reasonably homogeneous features.

In one embodiment, a geographically diverse region image such as a geographical image of United States of America can be divided into, for example, a plurality of polygonal image sections or regions (e.g., triangular sections, quadrilateral sections, pentagonal sections, hexagonal sections, etc., or any combination of two or more thereof) that have reasonably visually similar geographical features. Consider a geographically contiguous region. Also consider that when looking at this region from an elevated point (for example, of over 10,000 meters), the color-space of the region seems to be subjectively visually uniform (for example: a desert, a mountain range, or a tropical rain forest). In one embodiment, such a region can be considered as having reasonably similar features in the present context. Such a region can also be referred to as a relatively homogenous region in the present context. Histograms and LUTs can be computed for each of the plurality of polygonal image sections. One of the LUTs computed for a specific image section can be selected and applied to pixels that fall into that image section and having a higher bit-depth per band to transform each respective image section into a corresponding image section having a lower bit-depth per band that can be displayed or viewed on a display device. In other words, a specific look up table within the computed plurality of LUTs that is associated with the polygonal section to which a pixel with the higher bit depth per band belongs to can be selected and used to transform the pixel having the higher bit depth per band into a pixel having a lower bit-depth per band that can be displayed or viewed on a display device. A LUT generated by using the overall image histogram can be applied to pixels that fall in sections outside of any of the polygonal image sections having the features. However, the LUT generated by using the overall image histogram may not be needed if each pixel in the image falls within at least one specific polygonal section containing a feature. In other words, the LUT generated by using the overall image histogram may not be needed if substantially the entire image is divided into polygonal sections for which a specific LUT can be calculated.

For instance, for a given NUI image or R-set, the NUI can be divided into R distinct relatively homogenous regions. For each specific relatively homogenous region R_(s) (i.e., for each region having relatively homogeneous features), an entry E_(s) can be computed. The entry E_(s) is an abstract object that contains all statistical and polygonal data to associate E_(s) with a corresponding region R_(s).

In one embodiment, the polygonal sections can be quadrilateral sections. In one embodiment, an entry E_(s) of each quadrilateral section in the plurality of quadrilateral sections contains: (a) a single quadrilateral (quad) Q_(s). The quad Q_(s) is defined, for example, as a closely convex quadrilateral in pixel space to cover a region R_(s) that is reasonably homogenous with respect to its visual features F_(s). In one embodiment, the selection of the quadrilateral sections or regions having the homogenous features in the image can be accomplished manually by selecting manually each quadrilateral region or section in the image based on the perception of a user. In another embodiment, the selection of the quadrilateral sections can be accomplished automatically using methods or algorithms (implemented by hardware or software) that can detect various homogenous features and thus delimit the various homogenous features by quadrilaterals to define the quadrilateral sections.

In one embodiment, the entry E_(s) further contains (b) statistical data, for example in the form of a histogram H_(s) for all pixels within the quadrilateral Q_(s). The histogram H_(s) is computed and stored for each band in each pixel within the homogenous quad Q_(s) section. In one embodiment, pixels on the boundaries of each quadrilateral section Q_(s) can be included within the quadrilateral sections so as to eliminate, for example, sub-pixel errors or lost pixel issues. The relative homogeneity of each quadrilateral region or section is substantially not affected by the addition of the boundary pixels. The time t_(s) at which H_(s) was generated is also stored.

In one embodiment, the statistical data may further include, optionally, (c) a plurality of N computed LUTs, i.e., N LUTs L[N]_(s), where N is selected between zero and N_(max)(i.e., 0<N<N_(max)). N_(max) may be user-defined. The plurality N LUTs are generated and saved by applying methods that employ the histogram (e.g., histogram equalization, standard deviation, min-max stretching, linear contrast curves, gamma-based contrast curves, S-shaped contrast curves, Gaussian curves, etc.). The time t[N]_(s) at which each of LUTs L[N]_(s) were generated may be saved.

In one embodiment, the statistical data may also include, optionally, (d) a plurality of M custom LUTs L[M]_(s), where M is selected between zero and M_(max) (i.e., 0<M<M_(max)). M_(max) may be user defined. The time t[M]_(s) at which the LUTs were generated may be saved. A custom LUT is a LUT generated manually or derived empirically with manual intervention and control.

Hence, in one embodiment, an entry E_(s) is associated with each quad Q_(s) that covers a region R_(s) having homogenous features F_(s). Therefore, there are several entries E_(s) per NUI image as well as for each of its R-sets. A histogram H_(s) is computed for all pixels within a quad Q_(s) and stored at time t_(s). Optionally, a plurality of N-computed LUTs L[N]_(s) and a plurality of M-custom LUTs L[M]_(s) may be computed and stored at, respectively, time t[N]_(s) and time t[M]_(s). Therefore, E_(s) is an object whose members are Q_(s), H_(s), t_(s), and optionally L[N]_(s) and t[N]_(s), L[M]_(s) and t[M]_(s). Given a subjectively homogenous visual feature F_(s) that covers a geographic region R_(s), the corresponding entry E_(s) contains data that is sufficient to quantitatively define the geographic bounds of region R_(s) and convert all pixels that are within region R_(s) from a known bit-depth higher than 8-bits per band to 8-bits per band for visualizing these pixels on a 24-bit RGB display device. In one embodiment, E_(s) may be represented as E_(s)(Q_(s), H_(s), t_(s), L[N]_(s) at t[N]_(s), L[M]_(s) at t[M]_(s)).

In one embodiment, if the number of regions for a given image is relatively large, optionally, a bounding quad Q_(b) is generated for a group of T quads Q_(s)[T] (where T≧0) based on geographical proximity. In one embodiment, under a single bounding Q_(b), multiple quads Q_(s) may overlap. In addition, in one embodiment, multiple bounding quads Q_(b) may also overlap. FIG. 5 represents a hierarchical tree or relationship structure between the various quads, according to en embodiment of the present invention. The root of the tree or structure represents the entire image. As shown in FIG. 5, an image is divided into four regions R₁, R₂, R₃ and R₄, each defined by quads Q_(s1), Q_(s2), Q_(s3) and Q_(s4). The quads being defined, the image is further divided into two further quads Q_(b1) and Q_(b2). As shown in FIG. 5, Q_(b1) covers Q_(s2) completely and Q_(s1) and Q_(s3) partially. Also as shown in FIG. 5, Q_(b2) covers Q_(s4) completely and Q_(s1) and Q_(s3) partially.

In one embodiment, a hierarchical tree of overlapping bounding quads Q_(b) can be generated such that for a given bounding quad Q_(b):

(a) the siblings of the bounding quad Q_(b) are in the general proximity of the bounding quad Q_(b);

(b) the bounding quad Q_(b) wholly or partially envelopes the children Q_(s) of the bounding quad Q_(b); and

(c) the parent of the bounding quad Q_(b) wholly envelopes the bounding quad Q_(b).

In one embodiment, the root Q_(b(root)) of the hierarchical tree of overlapping plurality of bounding quads Q_(b) is defined as the bounding box of the NUI file, and the leaf nodes at the bottom of the hierarchical tree of overlapping plurality of bounding quads Q_(b) are all entries E_(s), each entry E_(s) being associated with a corresponding quad Q_(s).

For example, a single root Q_(b(root)) having a plurality of entries E_(s) (associated with quads Q_(s)) as children can be used. When a search is performed to find a quad Q_(s) containing a pixel P_(i) (see, FIG. 5), any well known tree-based search algorithm may be used. Consider the tree structure depicted in FIG. 5. If a pixel is in the region defined by Q_(s)2, for example, the search would look for the pixel from the root node Q_(b(root)). From the root node Q_(b(root)), the search would move to node Q_(b1) because the pixel is within the bounds of Q_(b1). Within bounding quad Q_(b1), the sections quads Q_(s)1, Q_(s)2 and Q_(s)3 are stored sorted in ascending or descending order of their geographical coordinates. By applying a known search algorithm through these regions, the pixel that is indeed within the quad Q_(s)2 can be found or located. Once this is known, the selected LUT from the corresponding entry E_(s)2 (associated with the quad Q_(s)2) is applied to convert P_(i) to the corresponding pixel in the lower bit-depth output image.

FIG. 6 depicts a structure for storing statistical data for the original image and each R-set, each divided into several regions, according to an embodiment of the present invention. This method of generating and storing feature-based statistics for a given image is applied to the original image and all its R-sets within the NUI file. The original image and each of its R-set is divided into several regions, each region R_(s) for the image or R-set containing an entry E_(s). The original image and each R-set also store the hierarchical tree to further assist in matching a pixel to a corresponding region quad Q_(s) and hence to the entry E_(s) containing Q_(s).

FIG. 7 is a flow chart depicting a method of pre-computing the statistical data used for converting a first image of higher bit depth to a second image of lower bit depth, according to an embodiment of the present invention. A NUI file containing an original image and a plurality of R-sets covering a large geographical area containing a variety of visual features, the visual features covering a wide gamut of color-space, are input, at S10. A first method or procedure (method A) is applied to the inputted NUI file, at S20. The method A includes dividing the image (original image) into visually homogenous regions, as described in the above paragraphs, at S22. Pixels within the homogenous regions may overlap. Pixels included in all regions cover or form the entire image. The method A further includes applying a sub-procedure or method B to each region R at S24. The sub-procedure (method B) will be described in detail in the following paragraphs. The method or procedure A further includes determining whether the number of regions within the image is a large number, and if the number of regions within the image is a large number, optionally computing a hierarchical tree to organize all regions into a tree-based structure, at S26. The method A further includes, storing the hierarchical tree-based structure into the file containing the image data, at S28.

The sub-procedure method B includes storing bounds of a region in a quad Q at S241. The sub-procedure B further includes computing histogram H_(s) of all bands for all pixels contained within the quad Q_(s), and record a time t_(s) at which it was computed, at S242. The sub-procedure method B further includes, optionally computing a number N of look-up tables or LUT_(s) (L[N]_(s)), each LUT containing a pre-computed pixel conversion table for all bands, to convert pixels at a higher-bit depth to a lower-bit depth, and record times at which the LUTs were generated, at S243. The sub-procedure method B may further include, optionally compute an additional number M of look-up tables or LUTs (L[M]_(s)), each containing a user-defined, custom pixel conversion table for all bands, to convert pixels at a higher-bit depth to a lower-bit depth, and record times at which the LUT_(s) were generated, at S244. The sub-procedure method B proceeds by creating an entry E_(s) containing Q_(s), H_(s), and t_(s), L[N]_(s) and t[N]_(s), and L[M]_(s) and t[M]_(s), at S245. The sub-procedure method B further stores the created entry Es into the file containing the image data.

FIG. 8 is a flow chart depicting the method of converting a first image to a second image that can be displayed for a given area of interest (AOI) or viewport in the NUI file, using feature-based statistics and associated data generated earlier, according to an embodiment of the present invention. The method includes, for a NUI file having a pixel depth per band greater than 8-bit per band, selecting a specific LUT from either L[N] or L[M] where N or M is known, at S30. The method further includes, reading pixels from an area of interest (AOI) or viewport from a NUI file, at S32. The areas of interest AOI or viewport can be for example a rectangular area within the bound quad Q_(b1), as shown in FIG. 5.

A viewport can be defined as a “window” within an image having a pixel width W and a pixel height H. The image can have a large pixel width (e.g., 1,000,000 pixels) and a large pixel height (e.g., 2,000,000 pixels). Such an image cannot be displayed on current display devices (e.g., a screen) in its entirety. In addition, such an image cannot be loaded in computer main memory in its entirety for any form of processing, given current available computer main memory capacity. Hence, an application is forced to load a smaller sub-section, window, area-of-interest or viewport of the very large image. This smaller sub-section, window or viewport has a pixel width VW and pixel height VH which typically range from several tens of pixels to several thousands of pixels (where W>VW and H>VH). For instance, within a 1,000,000 pixel wide and 2,000,000 pixel tall image, a sub-section or viewport starting from pixel coordinate (300000, 450000) within the image and having a pixel width of 1920 pixels and a pixel height of 1080 pixels can be selected. The viewport is wholly contained within the image. Therefore, if (0, 0) is the address of the top left pixel of the W×H image, and if (X, Y) is the address of the top left corner pixel of the VW×VH viewport within the image, then 0≦X<(X+VW)≦W and 0≦Y<(Y+VH)≦H.

The method further includes for each pixel P_(i), traversing Q_(b(root)) to find a Q_(s) that contains P_(i), at S34. If the pixel P_(i) is within multiple overlapping Q_(s), the Q_(s) having the most recent time t_(s) can be selected. The selected Q_(s) is associated to an entry E_(s). The entry E_(s) contains the type of LUT that the user selected. The method further includes, for each band in the pixel, determining its 8-bit value and populating the corresponding pixel in the converted image (i.e., converted AOI or viewport), at S36. In one embodiment, the selection of the Q_(s) for each pixel and the determination for each band of the pixel its 8-bit value is performed using parallel processing. However, as it can be appreciated, in another embodiment, this can also be performed using sequential processing. After each pixel is processed as described above, the converted image can then be displayed on a display device (e.g. a display screen), at S38.

The data that is generated using the above outlined procedure can be stored within a NUI file as a NUI Property. In one embodiment, the data can be formatted in Extensible Markup Language (XML) for readability. Alternatively, the data can be stored as a data structure in binary format.

For example, in the case of the continental United States of America (CONUS), if one were to generate multiple NUI image files to cover the entire continental United States, at a resolution of 1 meter per pixel, each NUI can range from few hundred gigabytes (GB) to several terabytes (TB). FIG. 1 depicts the geographical map of the continental United States of America (CONUS). An example of areas covered by the multiple NUI image files are shown as dotted rectangles on the geographical map. On the other hand, if the continental United States is divided into visually homogenous regions or sections (i.e., divided on regions having substantially homogenous features), quadrilateral shaped regions (quad regions) can be defined and shown on the geographical map of the continental United States of America. As shown in FIG. 1, the quadrilaterals can have various shapes including trapezoid shape, rectangular shape, or any other polygonal shape (e.g., pentagonal shape, hexagonal shape, or other more complex shape, etc.). As a result, each of the rectangles defining the NUI images can be divided into several quad regions. For example, as shown in FIG. 1, the NUI image of the area A₁ is divided into about 10 homogenous regions or sections R₁, R₂, . . . , R₁₀.

For each specific homogenous region R₁, R₂, . . . , R₁₀, a corresponding entry E₁, E₂, . . . , E₁₀) is available and stored in the NUI file. Any given pixel P_(i) (x, y) in a NUI file may fall within multiple regions and thus may be associated with a plurality of entries E_(s). For example, as shown in FIG. 1, a pixel may fall within homogenous section R₁ and within homogenous section R₃. In other words, the homogenous sections R₁ and R₃ overlap. As a result, the pixel that falls within both R₁ and R₃ may be associated with both entry E₁ and entry E₃. In this case, a simple rule may be applied to identify which entry to select (E₁ or E₃). For example, the entry E₃ can be selected if time t₃ is greater than time t₁.

A plurality of LUTs (N LUTs) can be generated by applying methods that employ a histogram H_(s) for all pixels within a given quad. The LUTs are generated using various methods including, histogram equalization, standard deviation, min-max stretching, linear contrast curves, gamma-based contrast curves, S-shaped contrast curves, Gaussian, etc. The N LUTs are then saved in the NUI file. For example, a plurality of LUTs can be generated for region R₁ by employing histogram H₁ within entry E₁. The plurality of LUTs generated for region R₁ are saved the NUI file that contains the region R₁.

A LUT determines the “look” of the image, after it is converted. The user may manually select what kind of LUT best describes the way in which the user wants to display the image prior to performing the conversion.

In the case where some pixels within a NUI file are replaced by other pixels, E_(s) entries for the regions containing the updated pixels also need to be updated accordingly. In another embodiment, the histogram update can be accomplished by re-computing the histograms H_(s) and the LUTs in the entry E_(s) for the updated pixels.

As stated in the above paragraphs, in one embodiment, the homogenous regions containing substantially homogenous features can be generated or determined manually based on the perception of a user or determined automatically using feature-detection methods or algorithms. In one embodiment, a database of regions that have substantially homogenous geographic features can be generated manually for the whole globe. The database can then be maintained or updated infrequently to take into account, for example, seasonal changes (spring, summer, fall and winter). Indeed, regions that were, for instance, greener during summer can be covered with snow during winter. As a result, regions that were part of one or more “greener quads” may be annexed to one or more “whiter quads.” In one embodiment, thermal and topological contour maps can be used to extend automation of the process of generating quads with homogenous geographic features. Desired regions having substantially homogenous features can be selected from the database and can be used to generate a list of histogram entries for a given NUI file.

FIG. 2A depicts a converted image from an unsigned 16 bits per band. The converted image contains a significant dynamic range of data, as shown on the top left area and on the bottom right area of the image. The image is converted using a conventional standard-deviation applied to the histogram of the entire image file. A LUT is derived and applied to convert the unsigned 16 bits per band image from the higher bit-depth per band (16 bits per band) to a lower bit-depth per band (8 bits per band) to display or view the converted image on a display device, the display device being configured to display an image having a pixel depth of 8 bits per band. This image shows what the output 8-bit depth image would look like without applying the method of converting an image having a higher-bit depth per band into an image having a lower bit depth per band described herein.

FIG. 2B depicts the same image converted from unsigned 16 bits per band to an unsigned 8 bits per band image according to an embodiment of the present invention. The regions have been shown in FIG. 2C. A histogram was computed for each of the regions. For each region, using the histogram, the corresponding unsigned 16 bit to unsigned 8 bit look-up table was generated using a conventional standard deviation method. The look-up table was then applied to the image within the corresponding region to convert from the original unsigned 16-bit image to a displayable unsigned 8-bit image.

The image shown in FIG. 2B looks “patchy” relative to the image shown in FIG. 2A which is relatively smooth and even. However, the “patchiness” of the image in FIG. 2B is due to the fact that the image in FIG. 2B is divided into polygonal regions (e.g. quadrilaterals) having substantially homogenous features so as to improve local detail by adjusting the dynamic range in each of the regions using statistical information computed for only those pixels that fall within those regions.

FIG. 2C depicts the polygonal layout (e.g., quadrilateral layout) used to divide the image in FIG. 2B into a plurality of regions or sections, according to an embodiment of the present invention. A plurality of LUTs are computed using histograms determined for each of the plurality of regions. As a result, the regions that define the “patchy” image display detail that is better than if a single LUT was derived from the histogram of the entire image, as is the case for the image shown in FIG. 2A.

FIG. 3A depicts an image portion obtained when zooming into the top-left corner region of the image shown in FIG. 2A. FIG. 3B depicts an image portion obtained when zooming into the top-left corner region of the image shown in FIG. 2B. The image portion shown in FIG. 3A lacks contrast and detail of the features in this image portion can not be discerned. On the other hand, the image portion shown in FIG. 3B and covering the same area covered by the image shown in FIG. 3A, is crisp and details of features in this image portion can be discerned.

Hence, the method of computing a histogram for each region in a plurality of regions allow for accommodating a greater dynamic range of color when converting pixel data from higher bit depths per band to lower bit depths per band (e.g., 8 bits per band) for display on a standard 24 bit RGB displays. By pre-computing and storing feature-based statistics or statistical data and histograms within images with the higher bit depths per band, conversion of an area or region of interest from a higher bit depth per band into lower bit depth per band can be directly performed (i.e., performed “on-the-fly”). This would eliminate the need to generate additional for-display-only 8 bits per band image conversion. By eliminating the need to generate additional for-display-only 8 bits per band image conversion, the use of storage disk space can be reduced or minimized.

In some embodiments, programs for performing methods in accordance with embodiments of the invention can be embodied as program products in a computer such as a personal computer or server or in a distributed computing environment comprising a plurality of computers. The computer may include, for example, a desktop computer, a laptop computer, a handheld computing device such as a PDA, etc. The computer program products may include a computer readable medium or storage medium or media having instructions stored thereon used to program a computer to perform the methods described above. Examples of suitable storage medium or media include any type of disk including floppy disks, optical disks, DVDs, CD ROMs, magnetic optical disks, RAMs, EPROMs, EEPROMs, magnetic or optical cards, hard disk, flash card (e.g., a USB flash card), PCMCIA memory card, smart card, or other media. Alternatively, a portion or the whole computer program product can be downloaded from a remote computer or server via a network such as the internet, an ATM network, a wide area network (WAN) or a local area network.

Stored on one or more of the computer readable media, the program may include software for controlling both the hardware of a general purpose or specialized computer or processor. The software also enables the computer or processor to interact with a user via output devices such as a graphical user interface, head mounted display (HMD), etc. The software may also include, but is not limited to, device drivers, operating systems and user applications.

Alternatively, instead or in addition to implementing the methods described above as computer program product(s) (e.g., as software products) embodied in a computer, the method described above can be implemented as hardware in which for example an application specific integrated circuit (ASIC) can be designed to implement the method or methods of the present invention.

Although the various steps of the method(s) are described in the above paragraphs as occurring in a certain order, the present application is not bound by the order in which the various steps occur. In fact, in alternative embodiments, the various steps can be executed in an order different from the order described above.

Although the invention has been described in detail for the purpose of illustration based on what is currently considered to be the most practical and preferred embodiments, it is to be understood that such detail is solely for that purpose and that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover modifications and equivalent arrangements that are within the spirit and scope of the appended claims. For example, it is to be understood that the present invention contemplates that, to the extent possible, one or more features of any embodiment can be combined with one or more features of any other embodiment.

Furthermore, since numerous modifications and changes will readily occur to those of skill in the art, it is not desired to limit the invention to the exact construction and operation described herein. Accordingly, all suitable modifications and equivalents should be considered as falling within the spirit and scope of the invention. 

What is claimed:
 1. A method for retrieving an area of interest of an image stored in a storage device, the image being associated with statistical data stored along with the image, the method being implemented by a computer system that includes one or more processors configured to execute computer program modules, the method comprising: selecting, by the one or more processors, the area of interest having user desired features, the stored image being divided into a plurality of polygonal image sections, each polygonal image section having substantially homogenous features, and at least two polygonal image sections in the plurality of polygonal image sections have different polygonal shapes defining different area sizes, each polygonal image section comprising a plurality of pixels and each pixel having a plurality of bands and a higher bit depth per band with a bit depth greater than 8 bits per band, the statistical data comprising statistical data of each band of each pixel within each polygonal image section, the statistical data being stored with the data of each corresponding polygonal image section in the storage device, the image being further divided into a plurality of quads such that each quad has boundaries that contains at least a portion of a first polygonal image section and at least a portion of a second polygonal image section; determining, by the one or more processors, in which quad in the plurality of quads each pixel of the area of interest is located; determining, by the one or more processors, for each quad, in which portion of the first polygonal image section or portion of the second polygonal image section each pixel of the area of interest is located; and upon locating each pixel with the higher bit depth per band within the area of interest, retrieving pixels of the area of interest from the storage device.
 2. The method according to claim 1, wherein the statistical data of each band of each pixel within each polygonal image section comprises histograms of each band of each pixel within each polygonal image section.
 3. The method according to claim 1, further comprising: converting, by the one or more processors, each pixel with the higher bit depth per band in the area of interest into a corresponding pixel of a lower bit depth per band using the statistical data of a polygonal image section in which each pixel is located to obtain an image corresponding to the area of interest having pixels of the lower bit depth per band.
 4. The method according to claim 3, further comprising displaying, by the one or more processors, the obtained image having pixels of the lower bit depth per band on a display device.
 5. The method according to claim 3, further comprising transmitting, by the one or more processors, the obtained image having pixels of the lower bit depth per band to a computer of the user in communication with the computer system upon request by the user.
 6. The method according to claim 3, wherein the obtained area of interest having the lower bit depth per band includes an image having a bit depth per band less than or equal to 8 bits per band.
 7. The method according to claim 1, wherein the image, the plurality of quads, and the plurality of polygonal image sections are organized in a form of a hierarchical tree structure such that a root node corresponding to the image has a plurality of branches corresponding to the plurality of quads and a branch node of the tree structure corresponding to a quad in the plurality of quads has two or more leaf nodes, each leaf node corresponding to a polygonal image section.
 8. The method according to claim 7, wherein determining, by the one or more processors, in which quad in the plurality of quads each pixel of the area of interest is located, and determining, by the one or more processors, for each quad, in which portion of the first polygonal image section or portion of the second polygonal image section each pixel of the area of interest is located comprises traversing the hierarchical tree structure from the root node corresponding to the image to the leaf nodes corresponding to the polygonal image sections.
 9. A computer program product comprising a non-transitory computer readable medium having instructions stored thereon when executed by a computer system performs the method recited in claim
 1. 10. A computer system of retrieving an area of interest of an image stored in a storage device, the image being associated with statistical data stored along with the image, the computer system comprising one or more processors configured to: select the area of interest in the stored image, the area of interest having user desired features, the stored image being divided into a plurality of polygonal image sections, each polygonal image section having substantially homogenous features, and at least two polygonal image sections in the plurality of polygonal image sections have different polygonal shapes defining different area sizes, each polygonal image section comprising a plurality of pixels and each pixel having a plurality of bands and a higher bit depth per band with a bit depth greater than 8 bits per band, the statistical data comprising statistical data of each band of each pixel within each polygonal section, the statistical data being stored with the data of each corresponding polygonal image section in the storage device, the image being further divided into a plurality of quads such that each quad has boundaries that contains at least a portion of a first polygonal image section and at least a portion of a second polygonal image section; determine in which quad in the plurality of quads each pixel of the area of interest is located; determine, for each quad, in which portion of the first polygonal image section or portion of the second polygonal image section each pixel of the area of interest is located; and upon locating each pixel with the higher bit depth per band within the area of interest, retrieve pixels of the area of interest from the storage device.
 11. The computer system according to claim 10, wherein the one or more processors are configured to convert each pixel with the higher bit depth per band within the area of interest into a corresponding pixel of a lower bit depth per band using the statistical data of a polygonal image section in which each pixel is located to obtain an area of interest having pixels of the lower bit depth per band.
 12. The computer system according to claim 11, wherein the one or more processors are configured to display or transmit an image corresponding to the area of interest with the pixels of the lower bit depth per band.
 13. A method for storing statistical data for pixels of an image within the image in a storage device, the pixels of the image having a bit depth per band greater than 8 bits per band, the method being implemented by a computer system that includes one or more processors configured to execute computer program modules, the method comprising: dividing, by the one or more processors, the image having pixels with the higher bit depth per band into a plurality of image sections, each image section in the plurality of image sections having a polygonal shape and including substantially homogenous features, at least two image sections in the plurality of image sections defining different area sizes, each image section comprising a plurality of pixels and each pixel having a plurality of bands and a bit depth greater than 8 bits per band; computing, by the one or more processors, statistical data for pixels within each polygonal image section, the statistical data comprising histograms of each band of each pixel within each polygonal section; and storing, by the one or more processors, the statistical data for pixels within each polygonal image section with the data of each corresponding polygonal image section in the storage device.
 14. The method according to claim 13, further comprising storing, by the one or more processors, data of each polygonal image section including storing coordinates of each polygonal image section in the storage device.
 15. The method according to claim 13, wherein the substantially homogenous features in each polygonal image section are determined based on a perception of a user or by a computer algorithm.
 16. The method according to claim 13, wherein the plurality of polygonal image sections comprise one or more of triangular image sections, quadrilateral image sections, pentagonal image sections, or hexagonal image sections, or any combination of two or more thereof.
 17. The method according to claim 13, further comprising further dividing, by the one or more processors, the image into a plurality of quads such that each quad has boundaries that contains at least a portion of a first polygonal image section and at least a portion of a second polygonal image section in the plurality of polygonal sections.
 18. The method according to claim 17, organizing the image, the plurality of quads, and the plurality of polygonal image sections in a form of a hierarchical tree structure such that a root node corresponding to the image has a plurality of branches corresponding to the plurality of quads and a branch node of the tree structure corresponding to a quad in the plurality of quads has two or more leaf nodes, each leaf node corresponding to a polygonal image section.
 19. The method according to claim 18, further comprising storing the hierarchical tree structure along with the image in the storage device.
 20. A computer program product comprising a non-transitory computer readable medium having instructions stored thereon when executed by a computer system performs the method recited in claim
 13. 21. A computer system for determining statistical data for pixels of an image within the image, the pixels of the image having a bit depth per band greater than 8 bits per band, the computer system comprising one or more processors configured to: divide the image having pixels with the higher bit depth per band into a plurality of image sections, each image section in the plurality of image sections having a polygonal shape and including substantially homogenous features, at least two image sections in the plurality of image sections defining different area sizes, each image section comprising a plurality of pixels and each pixel having a plurality of bands and a bit depth greater than 8 bits per band; compute statistical data for pixels within each polygonal image section, the statistical data comprising histograms of each band of each pixel within each polygonal section; and store the statistical data for pixels within each polygonal image section with the data of each corresponding polygonal image section in the storage device.
 22. The system according to claim 21, wherein the one or more processors are configured to store data of each polygonal image section including storing coordinates of each polygonal image section in the storage device.
 23. The system according to claim 21, wherein the substantially homogenous features in each polygonal image section are determined based on a perception of a user or by a computer algorithm.
 24. The system according to claim 21, wherein the plurality of polygonal image sections comprise one or more of triangular image sections, quadrilateral image sections, pentagonal image sections, or hexagonal image sections, or any combination of two or more thereof.
 25. The system according to claim 21, wherein the one or more processors are configured to further divide the image into a plurality of quads such that each quad has boundaries that contains at least a portion of a first polygonal image section and at least a portion of a second polygonal image section in the plurality of polygonal sections.
 26. The system according to claim 25, wherein the one or more processors are configured to: organize the image, the plurality of quads, and the plurality of polygonal image sections in a form of a hierarchical tree structure, such that a root node corresponding to the image has a plurality of branches corresponding to the plurality of quads and a branch node of the tree structure corresponding to a quad in the plurality of quads has two or more leaf nodes, each leaf node corresponding to a polygonal image section, and store the hierarchical tree structure along with the image in the storage device. 